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Abstract 



Taking the triangle areas as independent variables in the theory of 
Regge calculus can lead to ambiguities in the edge lengths, which can 

■ be interpreted as discontinuities in the metric. We construct solutions 

to area Regge calculus using a triangulated lattice and find that on 
a spacelike hypersurface no such discontinuity can arise. On a null 
hypersurface however, we can have such a situation and the resulting 

■ metric can be interpreted as a so-called refractive wave. 

o 

O ' 1 Introduction 

The Ponzano- Regge ^ model of quantum gravity in three dimensions begins 
with a labelling of the 1-simplices in a simplicial complex with irreducible 
i , representations of SU(2). Using these labellings a state sum is constructed, 

£jpj| which is the discrete analogue of the partition function of three dimensional 

quantum gravity. The sum over labellings of the 1-simplices in the state 
sum corresponds to an integration over edge lengths in the path integral. 
The action appearing in the exponential of the path integral is the Regge 
action |2j, defined in terms of these edge lengths. 

In four dimensions the labelling of the 1-simplices does not seem to be 
sufficient j^j- For example in the Barrett-Crane model [3] we label the 2- 
simplices with balanced irreducible representations of SU(2) x SU{2). The 
resulting state sum then takes the form of the path integral for the four 
dimensional Regge action, yet now the areas of the triangles appear to be 
taking the role of the independent variables. 

This is then our motivation to study Regge calculus with triangle ar- 
eas, rather than edge lengths, as the independent variables. This idea was 
originally suggested by Rovelli 5 j in connection with loop quantum gravity. 

*e-mail address: C.J.Wainwright@damtp.cam.ac.uk 
^ e-mail address: R.M.Williams@damtp.cam.ac.uk 
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In Section 2 we will introduce the theory of area Regge calculus and 
discuss some if its potential consequences such as discontinuous metrics. In 
Section 3 we will construct full solutions of area Regge calculus to investigate 
the possibility of discontinuities along spacelike and null hypersurfaces. We 
then, in Section 4, need to introduce the theory of tensor distributions which 
allows us, in Section 5, to draw some interesting comparisons between area 
Regge calculus and the theory of gravity. 

2 Area Regge Calculus 

In area Regge calculus we begin with the Regge action on a 4-dimensional 
simplicial complex |2j, 



However, instead of taking the lengths of the 1-simplices to be the inde- 
pendent variables we consider the areas of the 2-simplices, i.e. the Aa, to 
be independent. Variation of the action with respect to the A a results in 
the field equations e/± = 0, which states that the deficit angle at each tri- 
angle is zero. If this were standard Regge calculus, a zero deficit condition 
would imply the space were locally flat. However, the situation here is less 
straightforward, as explained in j^j, which we shall now summarize. 

Consider a single 4-simplex, which contains 10 edges and 10 triangles. 
The edge lengths uniquely specify the geometry of the simplex and hence 
determine the areas. However the reverse is not true; two 4-simplices can be 
constructed having the same triangle areas, yet different edge lengths. This 
problem can be overcome if we consider our 4-simplices to be in some sense 
close to a regular 4-simplex, analogous to choosing the principal value of a 
multi- valued function such as sin -1 x. 

Now suppose we have two 4-simplices meeting at a tetrahedron. We 
now have a total of 16 triangles and 14 edges, hence a problem of over- 
determinism. Also note the shared tetrahedron has 4 triangles and 6 edges, 
hence an under-determined geometry. We can thus envisage a situation 
where calculating the edge lengths of the tetrahedron from the areas of one 
4-simplex results in a different answer to if they were calculated from the 
other 4-simplex. We could interpret such a situation as a discontinuity of 
the metric. 

We shall proceed by investigating under what circumstances such dis- 
continuities can arise in solutions to area Regge calculus. 




(2.1) 



A 
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3 Hypersurface Discontinuities 

We have seen in Section 2 how discontinuities might arise in area Regge 
calculus; however this was only for two adjoining simplices. We would like 
to be able to construct a full solution to area Regge calculus containing such 
discontinuities. 

The simplest situation would be to restrict the discontinuity to some 
hypersurface S, separating M into M^. We triangulate M such that we 
also triangulate S. Taking the areas of the triangles as the independent 
variables, we have the field equations €a = 0, i.e. zero deficit angles. 

To restrict the discontinuity to £ we demand the areas are chosen to 
give well defined edge lengths in respectively. Combined with the zero 
deficit condition we see that both M are flat. 

We simplify the situation further by assuming £ is embedded in 
with no extrinsic curvature; this takes care of the zero deficit condition by 
virtue of the fact that the interior dihedral angles at those triangles in £ 
will be 7r in both M^. We work with a triangulated lattice, each lattice site 
being identical to its neighbours except across £. 

The procedure we will take is as follows. 

• Begin with a flat 3-lattice in dM~ . 

• Transform the coordinates of the lattice points (thus keeping the lattice 
flat) in such a way that the areas of the triangles are unchanged. This 
gives the lattice on dM + . 

• Extend the 3-lattices on dM^ to 4-lattices on and since the tri- 
angle areas agree we can identify the lattices on 3M ± giving the full 
4-dimensional solution. 

We can look at the cases of £ spacelike and £ null separately. 
3.1 S Spacelike 

We work in coordinates (i, x, y, z) with £ located at t = 0. Our flat 3-lattice 
in dM~ is given by coordinates in (x,y,z). The metric in this space is 
ds 2 = dx 2 + dy 2 + dz 2 . We label the edges of the triangulated 3-lattice 
from 1 to 7 as shown in Figure I3~T1 consistent with the binary notation of 
Since the hypersurface is flat the coordinates of edges 3, 5, 6 and 7 
are determined by those of edges 1, 2 and 4. 

Let Xi be the vector along the i'th edge, then a transformation of the 
lattice is determined by a transformation of X\, X2 and X4. This induces a 
linear transformation on R 3 , which we write as a matrix T £ GL(3). 
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Figure 3.1: Triangulated lattice edges 



Now since each lattice site is identical we only have 6 independent trian- 
gle areas 1 , these are A±z, A15, A\j, A2Q, A27 and A37, where Aij is the area 
of the triangle which has i and j as edges. 

Under the transformation Xj — > x\ = Txi we require that these six 
areas remain invariant. We can decompose T in order to analyse which 
transformations will satisfy this condition. 

We proceed by first noting we can always split T as T = OP, where 
O € 0(3) and P is some positive definite symmetric matrix. Since P is 
symmetric we can further decompose it as P = QDQ T , where Q 6 0(3) 
and D is diagonal with positive entries (since P was positive definite). Thus 
we have decomposed T as 

T = OQDQ T . (3.1) 

Since both O and Q are orthogonal they will not transform the triangle 
areas. Thus we only require that D does not transform the triangle areas. 
So, let D = diag(a, b, c) with a,b,c> 0, and consider a triangle based at the 
origin, with the other two vertices at (Xa, Ya, Za) and (Xb, Yb, Zb). This 
triangle has an area A given by 



Y A Z a 
Y b Z b 



+ 



Za Xa 
Zb Xb 



+ 



X A Y A 
X B Y b 



(3.2) 



1 Each lattice site has 12 triangles associated to it, however in the flat case each of these 
triangles has another which is always congruent to it, for example ^13 = A23, An = Aq7, 
etc. 
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and under a transformation by D the area changes as 



1 

4 





Y A 


Z a 


2 

i 2 2 
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Z a 


X a 


2 

+ aV 
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Y b 


Z b 


Z B 


X b 


X b 


Y b 





(3-3) 

Now, we require A' = A. Hence, if we let a = b 2 c 2 — 1, j3 = c 2 a 2 — 1 and 
7 = a 2 b 2 — 1 we must have 



a 



Y A Z a 
Y b Z b 



+ 13 



Za Xa 
Zb Xb 



+ 7 



X A Y a 
X b Y b 



0. 



(3.4) 



This condition has to be satisfied for all 6 triangles, giving 6 linear equations 
in a, P and 7. Letting v T = (a,(3,~f) we can write these in matrix form, 



Gv = 0, (3.5) 

where G is a 6 x 3 matrix, each row consisting of the squared determinants 
appearing in equation ()3.4|) for each of the six triangles. 

We argue in appendix IA.1I that the matrix G will have rank 3. Hence 
the only solution to equation ()3.5|) will be v = 0. This in turn implies that 
a = b = c = 1 and thus only orthogonal transformations leave the 6 triangle 
areas independent. 



3.2 E Null 

Now we wish to perform a similar analysis for the case when T, is null. We 
work in double null coordinates (u, v, y, z) with E located at u = 0. Our flat 
three lattice is given by coordinates in (v,y,z). The metric in this space is 
simply ds 2 = dy 2 + dz 2 , since v is a null direction. 

Again we transform the lattice by changing the three linearly indepen- 
dent vectors defining it; this induces a GL(3) transformation on the (v , y, z) 
space, which is given by the matrix T. We proceed to decompose this trans- 
formation, as we did in the spacelike case, but in a different fashion. Note 
now that areas of triangles do not depend on the v coordinates of their 
vertices, hence we begin by writing T as 



fa b T \ _ ( a- b T E- x c b T E- 1 
V c E )~\ 1 

where E € GL{2) and b, c are 2-component vectors. We can decompose 
E as we did for T in the spacelike case, giving E = OP = OQDQ T , with 
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0,Q € 0(2) and D diagonal with positive entries. Hence T now takes the 
form 



* 





OQ 



which can be rearranged to the form 



OQ 



1 
d 





b 



o 

D 





Q T 



(3.7) 





Q T 



(3.8) 



where d = (OQ) T c, and elements denoted by an asterisk (*) will not be 
relevant to our calculations. 

We consider a triangle based at the origin with its other vertices at 
(Va, Ya, Z a) and (Vb, Yb, Zb)- The area of the triangle is given by 



A 



Y A Z a 
Y B Z B 



(3.9) 



which does not depend on the v coordinates of the vertices, as stated earlier. 
We see that the first and third matrices in the decomposition of equation 
(13.811 will not alter the areas of the triangle, so all we require is that the 
second matrix does not either. 

Let d T = (c, d) and D = diag(a, b) with a, b > 0, then what we require is 
that the triangle areas remain invariant under y' = ay + cv and z' = bz + dv. 
For our triangle's coordinates this gives the condition 



(ab-1) 



Y A 


Z a 


-be 


Z A 


V a 


— ad 


Ya 


Y a 


Y b 


Z b 


Z b 


Vb 


V b 


Y b 



0. (3.10) 



Again we demand that the 6 triangle areas, j4i3, A15, An, A26, A27 and 
A37, remain invariant under the transformation. Thus we get 6 equations in 
(ab — 1, —be, —ad), which we can write as 



our 3 unknowns v T 



Gv = 0, 



(3.11) 



where G is the 6x3 matrix whose six rows contain the determinants ap- 
pearing in eauation ()3.1U|) for each of the six triangles. We show in Appendix 
IA. 21 that the matrix G has rank 3, and hence the only solution to equation 
(|3.11|) is v = 0. Thus we have ab = 1 and c = d = 0, which leads to the 
most general form of T leaving all six triangle areas invariant as 



* * 
OQ 



1 
D 



1 

Q T 



(3.12) 
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where again * denotes elements which will not be relevant to us, and the 
matrix D takes the form 

D= ( I °, ) • (3.13) 

with a > 0. 

3.3 Metric discontinuities 

We can interpret these solutions to area Regge calculus as defining a metric 
on the space M with a possible discontinuity at E. We began with a lattice 
in E, defined on a metric which we shall call g. The edge lengths of the lattice 
were defined in terms of this metric by I 2 = x T gx. We then transformed the 
lattice with a linear transformation, which sent x — > Tx. Thus the lengths 
of the edges were transformed as I 2 — > I' 2 = (Tx) T g(Tx). Equivalently, we 
can think of this as defining a new metric on E given by 

9' = T T gT. (3.14) 

Edge lengths are then calculated using the old coordinates I' 2 = x T g'x. 

For the case of E spacelike, our initial metric was 5=13, and we found 
that our transformation T had to be orthogonal. Thus, the new metric is 
given by g' = T T t^T =13. So we see that in this case there is no difference 
in the metrics and the edge lengths will then be well defined across E. 

For the case of E null, our initial metric was of the form 

9=(l I), (3-15) 



and we found that T must take the form given in equation (|3.12|) . Using 
this we can calculate the form of the new metric, which turns out to be 

9 >= ( n nfflnT ) > ( 3 - 16 ) 



QD'Q 1 

where D takes the form given in equation (J3.13|) . More explicitly, letting Q 
be a rotation through angle 6, we have 




g'=\ a 2 cos 2 6 + aT 2 sin 2 9 (a 2 - a~ 2 ) cos 9 sin 9 ). (3.17) 
(a 2 — a~ 2 ) cos 9 sin 9 a 2 sin 2 9 + a -2 cos 2 9 

This solution is precisely the metric of a refractive wave spacetime jS]. 
Note that if we take a to be close to unity, i.e. a = 1 + e/2, then to first 
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order the difference in metrics on the hypersurface will take the form 



/ \ 

g'- g = ei cos 26 sin 26 , (3.18) 
\ sin26> -cos 26 J 

and we can clearly see that the angle 6 can be interpreted as the polarization 
of this wave, and the value a as its strength. 

4 Tensor Distributions 

In order to handle discontinuous quantities properly we need to use distri- 
butions. The standard idea of a distributions can be extended to tensor 
valued distributions on arbitrary manifolds. There have been many differ- 
ent approaches to this idea lll[ I12| ; we follow here Geroch and Traschen 

us- 

Just as we can define distributions on R as linear maps from a space of 
test functions to R, so we define tensor distributions on a manifold M as 
linear maps from a space of test fields to R. The test fields in this case are 
tensor densities of weight +1 with compact support, and we write the linear 
map as 

which we need to satisfy certain continuity conditions, which are discussed 
more in |12j . 

Every smooth tensor field S£ "' on M defines a distribution by 

(st :,€::.) = J W^t.dv. (4.2) 

The fact that we use densities, instead of just tensors as the test fields means 
we do not need to introduce a volume element for this integral, as is done 

in manu. 

We can define the derivative of distributions, which is consistent with 
that for smooth tensors, by 

iiirt;. r;;;;;.n, <:/';-. (4.3) 

If we introduce a smooth connection on M then we have the corresponding 
result 

(v c f b a : ;:, Pt.) = -(7 6 a ::- : v,#t;;.)- (4.4) 
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We can use the derivative of distributions to define the weak derivative 
of locally integrable tensors 2 . A locally integrable tensor field T^" defines 
a distribution TP'", the same way as for a smooth tensor field in equation 
(|4.2|) . The weak derivative of T b °"" is - if it exists - a locally integrable tensor 
field R a cb such that R%- = V c f b a -. 

The curvature tensor is nonlinear in the metric, and since multiplication 
of distributions does not usually make sense, we need to investigate which 
distributional metrics lead to a well-defined curvature tensor. 

Choosing an arbitrary smooth connection V c , with curvature tensor 
P a bcd, we can write the curvature tensor for a smooth metric g a b as 

R a bcd = P a bcd + 2Q a e [cQ e d]b + 2V[cQ a d]b, (4-5) 

where 

Q\c = \g ae (2V(b9c)e - V e5 fe c ) • (4.6) 

Equation ()4.5|) continues to make sense if we relax the smoothness conditions 
on g a b. Geroch and Traschen show that a distributional curvature tensor can 
still be defined for a so-called regular metric. A metric is said to be regular 
if (i) g a i, and g ab exist everywhere and are locally bounded, and (ii) the weak 
first derivative of g a b exists and is locally square integrable. 

Now consider a discontinuous metric. We split the manifold as M = 
M + U M~ with M + n M~ = S, where S is a codimension 1 smooth hy- 
persurface. The metric is then given by smooth tensor fields g^ b defined on 
M respectively. This metric clearly satisfies condition (i); to calculate the 
weak first derivative we note that this metric defines the distribution 

g ab = e+ 5 + + e- g - b , (4.7) 

where are the step functions on M^ 1 . The first derivative of this is 

VcQab = e + V c <7+ + @-V c g- b + Snjgab], (4.8) 

where 5 is the delta function on the boundary dM~ , n c is the normal to 
this boundary, and [g ab \ = g^ldM- ~ 9ab\dM~- Since there is no locally 
integrable tensor field which gives rise to the delta function we see that if 
{dab} 7^ then the weak derivative of g a b does not exist, hence condition (ii) 
is not met, and the metric is not regular. 

It is worth pointing out that distributional valued metrics appear quite 
frequently in the literature, in particular with regard to impulsive gravi- 
tational waves. Occasionally it is more convenient to write the metric for 

2 A tensor field T£'[' is locally integrable if T^'"(f> h a y is Lebesgue measurable and its 
integral converges for all test fields <j> a ;;, . T^ "' is locally square integrable if T^ 'T^ y is 
locally integrable 
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such solutions in a coordinate system which is not continuous. When this 
is done the metric may for instance contain delta functions, making it un- 
bounded and hence not regular. However this is not a problem because 
the coordinate system was considered discontinuous, and a transformation 
of the coordinates back to a continuous coordinate system yields a regular 
metric. For a nice discussion of this see Here, we are assuming the 

coordinate system to be continuous, and as such the discontinuities cannot 
be transformed away with such discontinuous coordinate transformations. 

Another problem with a discontinuous metric is the ambiguity of parallel 
transport. This is due to the fact there is no well defined metric connection. 
The standard formula for the metric connection will fail due to ill-defined 
multiplication of distributions, so the best we can do is a so-called regularly 
discontinuous connection Here the connection components are contin- 
uous, and metric compatible, in yet not across the boundary, which 
results in the following rule for differentiation of the metric, 

V c 0a& = Sncfgabl (4.9) 

Despite the lack of a natural rule for parallel transport, we can still define 
metric geodesies. Since we are just extremizing the path length this problem 
is well defined. The result is familiar from optics, where we can describe the 
path of light through a medium with a metric representing the time of travel 
along the path. Where the refractive index of a medium changes abruptly 
we have a discontinuity of this metric, and we know that the path of light 
is refracted at this point. Exactly the same occurs here, with the geodesic 
refracted at the boundary by the rule 



g + (X + ,T) = g-(X-,T), VTgT p S, (4.10) 
g+(X + ,X + ) = g-(X-,X~), (4.11) 

where X^ is the tangent to the geodesic on either side of the boundary. 

It is interesting then that spaces which are solutions to area Regge cal- 
culus may not have a well defined notion of parallel transport. However it 
seems we would be able to define geodesies in the space. 

5 Comparisons with General Relativity 

We have shown in Section 4 that a discontinuity in the metric does not make 
sense in general relativity, since the curvature tensors cannot be defined. 
However, it is interesting to analyse the linearized Einstein equations, since 
these allow distributional solutions of all kinds. This will provide us with 
an interesting comparison with our area Regge calculus solutions. 
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We shall work with the perturbation about flat space g ab = f]ab + £h ab 
and write the linearized equations as 

Gf f cd h ab , cd = 0, (5.1) 

where Gf/ d = r] ac 5 b e 5 d f + rj ac S b f d d - if d b a e b) - r] ef r] ac r] bd is a constant tensor 
and h a b = h ab — \r} ab h is the trace reversed perturbation. 

Considering h ab as a tensor distribution (from now on we will drop hats), 
we can find solutions to equation (|5.1j) by its action on an arbitrary test 
tensor density (fi e f . 



(Gf f cd h abM ,^) = (5.2) 

^Gf f cd (h ab ,cf> e ^ cd ) = 0. (5.3) 

Substituting in the form of our desired discontinuous solution 

hab = <5> + Pab + @~q a b, (5-4) 

we obtain the integral 

Gf f cd [ P ab cp ef ,c d dV + (g-term) = 0. (5.5) 
Jm+ 

Integrating by parts twice gives 



/-iabcd 



+(g-term) = 0, 



Pah,cd4> ef dV + / {p ab n c (/) ef \ d -p ab c n d cf) ef )dS 

W+ JdM+ 

(5.6) 

where n a is the normal to the hypersurface S. However we should be careful, 
since the metric is discontinuous, the normal vector will also be discontinuous 
across the hypersurface. In fact we can define n\ = n a + em± and then solve 
for m|. Fortunately, as the discontinuity in the normal direction is 0{e) it 
will not come into our calculations, but will be important later. 

Now <fi e f was arbitrary, hence choosing supp((/> e ^) C M + only the first 
integral contributes, hence we require the linearized Einsteins equations to 
be satisfied in M + , i.e. 

Gf } cd P abM = 0. (5.7) 

Similarly q ab satisfies the linearized Einstein equations in M~ . 

The more interesting result comes from the surface integral. To deal 
with this we must split the derivative into parallel and perpendicular parts, 
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d c = dc +n c d . Integrating the parallel derivative components by parts we 
obtain 



/-yabcd 



'dM+ L 



{Vabncn d )d^(f S - {Pab,c n d + d\{p ab n c ))4) eS dS + (g-term) = 



(5.8) 

Note that we can choose (p e f |s and d^(\f^ |s independently hence we 
obtain two junction conditions satisfied on the surface S, given by 

Gf f cd uJ ab n c n d = (5.9) 
Gf f cd ((u ab n c ) 4 + Lj abtC n d -d ± (U ab n c )n d ) = 0, (5.10) 

where we have defined U ab = p ab |s — q ab |s and derivatives as uj abfi = p ab c |s 
— g a bcls- To interpret these conditions we look at two important cases. 

5.1 E Spacelike 

Junction conditions at spacelike surfaces in general relativity are well known 
and were given by Israel |14j : we show here that the above method repro- 
duces these results. 

We work in coordinates (x°, x 1 , x 2 , x 3 ) with rj ab = diag(l, — 1, — 1, — 1) 
and take S to be given by x° = 0, hence n c = (1,0,0,0). Equation (|5.9j) 
can be split in (0, 0),(0, i) and components, where ... = 1,2,3. 

The (0,0) and (0, i) components are automatically satisfied, and the 
components give the condition 

Wij = 0, (5.11) 

hence the induced metrics on S must agree. Note that there is no restriction 
on the (0, 0) and (0, i) components of the discontinuity. This is related to the 
discontinuity in the normal to the hypersurface, mentioned earlier. A direct 
calculation reveals that the normal takes the form n 1 ^ = (1 — ^e/ij^, e/i^) 
on each side of the hypersurface. Thus the discontinuity in the normal 
direction is equivalent to the discontinuity in the (0, 0) and (0, i) components 
of the metric. Note also that this discontinuity could be transformed away 
with a continuous (yet not C 1 ) coordinate transformation, though such a 
transformation will change the differentiable structure of the manifold. 

Equation ([5.10ft can also be split, with the components giving the 
only condition, which we can write as 

Wij,0 — ^OiJ — WQj,i — Vij ( u k,0 ~ 2w fc Q,fc) = 0. (5.12) 
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which is simply the linearized form of the Israel matching condition in the 
vacuum, 



\Kij - rHjK] = 0. (5.13) 

Hence the extrinsic curvatures, Kij, of S must agree from both sides. 

This result is then completely analogous to that in area Regge calculus 
where we found that we could not have a discontinuity in the metric, and also 
the extrinsic curvature had to match due to the zero deficit angle condition. 



5.2 E Null 

Junction conditions for null hypersurfaces have been studied in [151 116) : 
however all assume the induced metrics on the null hypersurface agree. Fol- 
lowing on from [5], and since we are only considering the linearized case, we 
do not make such an assumption, and see what the field equations demand. 

We work in coordinates (it, v, y, z) with metric ds 2 = 2du dv — dy 2 — dz 2 . 
Let £ be given by u = hence n c = (1,0,0,0), but note d 1 - = Jjj- = l c d c , 
where l c = (0, 1,0,0). Despite this the above integrations by parts can still 
be done, giving the same equations (|5.91 15.10")) . 

Let indices a, (3, ... represent the spatial (y,z) directions. Solving equa- 
tion ([5.9)1 . the only non-trivial solutions are the (u,u), (u,a) and (a, (3) 
components, which can be combined into the following one equation 

ZJ ab n b = 0. (5.14) 

This is exactly the condition obtained by Barrett which states that the 
areas of spacelike 2-surfaces is continuous across E. This is therefore the 
linearized metric of a refractive wave spacetime. Solving for w a b, we can see 
explicitly the form of the metric discontinuity. 



Uab 



( 






^uy 




\ 




^uv 















^uy 





LOyy 






V 









— LOyy 





(5.15) 



Note again that the LO ua terms can be transformed away with a continuous 
transformation, yet the LO a p terms cannot. 

This result is also completely analogous to that obtained in area Regge 
calculus, where we found a discontinuity which, in the linearized situation 
of Equation (|3.18j) , is precisely of this form. 

Solving equation (|5.1U|) and using equation 1)5.14(1 to simplify the results 
we obtain four conditions, which we write in the following suggestive forms 
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V (Uau,0 + ^Pu,a ~ Ua/3,u) = 0, (5.16) 

^ a,j3 + (UvattU ~ ^vu,a ~ ^ua,v) = 0, (5-17) 

^uv,v ^vu,v ~ ^vv,u = 0, (5.18) 

Wva,/3 + U v p ia — U a p,v = 0. (5.19) 



Now, for a null hypersurface we can define various fundamental forms, 
see 1 15] for more details. These are the internal second fundamental form 
XaPi the- external second fundamental form ip a p, the normal fundamental 
form 7] a , and a quantity u. In this case these are given by 



Xaf3 = Van/3, (5.20) 

V>a/3 = VaZ/j, (5.21) 

n a = VJ V , (5.22) 

u = V v l v , (5.23) 



which we can calculate for this linearized situation, and hence rewrite equa- 
tions 1)5. 16115.19]) in terms of these quantities 



0, (5.24) 

0, (5.25) 

0, (5.26) 

0. (5.27) 

The quantity T a ^ in the second of these equations does not appear to be 
interpretable in terms of our fundamental forms. 

In the case where the induced metric is assumed to be continuous these 
conditions reduce to those obtained by Clarke and Dray ]15j . namely that 
fn^ipafs} = frjal = [w] = 0. However, here we have the more general case 
allowing for the discontinuous metrics present in the linearized thoery. 

6 Conclusions and Outlook 

We have shown how to construct non-trivial solutions to area Regge calculus. 
In particular we have shown the existence of solutions with a discontinuity 
in the metric along a null hypersurface. These solutions are precisely the 
refractive wave spacetimes introduced by Barrett (Hj- 

The relation between our results for area Regge calculus and those of 
the linearized Einstein equations is also apparent. At spacelike hypersurfaces 
both theories predict a matching of the induced metrics, and of the extrinsic 



PV + 277J = 
M = 

had = 
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curvature. At null hypersurfaces, both theories predict the discontinuous 
metric of a refractive wave spacetime. We would also like to be able to 
compare the curvature matching conditions in this case. 

This comparison with linearized gravity fits in nicely with the ideas of 
Rocek and Williams They calculate the weak field expansion of area 
Regge calculus and show the dynamics of the theory to be equivalent to 
that of edge length variable Regge calculus. Since at the perturbative level, 
in the long wavelength limit, Regge calculus is equivalent to Einstein's theory 
[Z1IHI, so also will area Regge calculus. 

Despite this perturbative/dynamical relation, area Regge calculus is 
clearly not so simply associated with the full theory of general relativity. 
The refractive wave solution of area Regge calculus is not a solution of the 
full non-linear theory of relativity, since for such a metric the curvature 
tensor will not be well defined. 

The next step in the study of area Regge calculus would be to look for 
more complex and interesting solutions. It appears to be straightforward 
to combine refractive wave solutions so long as they are moving parallel or 
antiparallel, the result being a spacetime with multiple discontinuities. In 
double null coordinates (u, v, y, z) the discontinuities are along surfaces of 
constant u, and surfaces of constant v. These discontinuities will thus inter- 
sect, and we can check that the field equations are still satisfied where this 
occurs. We can also consider geodesies which start parallel, but travel ei- 
ther side of such an intersection to a new region of space. Here the geodesies 
will no longer be necessarily parallel, which reveals the non trivial nature of 
these solutions. 

More general solutions to area Regge calculus can presumably be found. 
Whether we can relate the properties of these solutions, such as geodesic 
behaviour, to properties of solutions to general relativity remains to be seen. 

The ultimate goal will be to relate area Regge calculus back to the quan- 
tum theories of gravity from which it was born, in order to understand to a 
greater extent what these theories are telling us about the nature of space, 
time and gravitation. 
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A Appendix 

A.l Rank of G 

We show here the matrix G of Equation (|3.5|) has rank 3. First we can define 
a vector sab by 



(A.l) 



°AB 

Then G takes the form 
G T 



( 


Y A 


Z a 


2 


z A 


X a 


2 


X A 


Y a 






Y b 


Z b 




Z B 


X b 


5 


X b 


Y b 





(Sl3, Si 5 , S17, S 26 , «27 , S37) 



(A.2) 



Now we know the coordinates of vertices 3, 5, 6 and 7 are given in terms 
of those of vertices 1, 2 and 4. So, using the properties of determinants 
we have various relations to simplify the matrix. We can then use column 
operations on the matrix G T to put it in the form 



G T = (S24,S41,S12, *,*,*) • 



(A.3) 



Without loss of generality we choose the coordinates of our lattice basis 
vector as x± = (a,b,c), xT^ = (0, d, e) and xj = (0,0,/) where a,d,f ^ 0. 
The matrix G will then take the form 



G = [ uu , (A.4) 



thus, since a,d, f 7^ 0, this matrix will have rank 3. 
A.2 Rank of G 

We show here that the matrix G of Equation (|3.11|) has rank 3. First define 
the vector tab by 



/ 


d 2 f 








\ 




* 











* 


* 


a 2 d 2 






* 


* 


* 






* 


* 


* 




V 


* 


* 


* 


/ 



T 
r AB 





Y A 


Z a 




Z A 


Va 




V A 


Y a 


( 


Y b 


Z b 




Z b 


V b 


5 


V b 


Y b 



Then we can write G in the form 



G T = (ri 3 , r 15 ,r 17 , r 26 , r 27 , r 37 ) . 



(A.5) 



(A.6) 



Now we know that the coordinates of the vertices 3, 5, 6 and 7 can be 
given in terms of those of vertices 1, 2 and 4. So, using the properties of 
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determinants we have for example 7*13 = v\\ + T\ 2 = r\ 2 , and likewise for 
the other entries of G. Thus we can use column operations to put G T in the 
form 

G T = (r24,r 4 i,ri2, 0,0,0). (A.7) 
The first three rows of G are then just the cofactor matrix of 

/ V 1 Y l Z l \ 
S = \ V 2 Y 2 Z 2 , (A.8) 

V v A y 4 z 4 J 

and since the three vectors along the edges 1, 2 and 4 are linearly inde- 
pendent we have det 5^0. Hence, the cofactor matrix of S will also be 
invertible, which implies that G will have rank 3. 
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